Ground states and magnetization process for an triangular lattice array of magnetic 

dots with perpendicular anisotropy 
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USA. 



We analyzed the ground state of the array of magnetic particles (magnetic dots) which form a 
two-dimensional triangular lattice, and magnetic moment of which is perpendicular to the plane of 
the lattice, in the presence of external magnetic field. In the small fields long range dipole-dipole 
interaction leads to the specific antiferromagnetic order, where two out of six nearest neighbors of 
the particle have the same direction of magnetization moment and four - the opposite one. It is 
shown that magnetization process in such array of particles as opposed to the rectangular lattices 
results from the formation of the magnetized topological defects (dislocations) in the shape of the 
domain walls. 

PACS numbers: 75.10.Hk, 75.50.Tt, 75.30.Kz 
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I. INTRODUCTION. 

Magnetic ordering is usually attributed with exchange 
interaction of atomic spins, leading to rather simple 
magnetically ordered states. 1,2 Long-range magnetic 
dipole interaction usually produces smooth non- 
uniformity (domain structures of different kind) above 
this simple exchange structured— Nevertheless, systems 
of magnetic moments with pure dipolar interaction, 
so-called dipolar magnets, have been theoretically 
studied for more than sixty years^ and many physical 
properties, lacking in the spin-exchanged systems, are 
known for those models. Note first the presence of 
a non-unique ground state with nontrivial continuous 
degeneracy for quite simple bipartite lattices, like three- 
dimensional cubic lattice,— and for two-dimensional 
square lattice^— as well as specific phase transitions 
induced by external magnetic fieldj^— Magnon spectra 
for dipolar magnets demonstrate non-analytic behavior 
either for small wavevectorsj^— or at some symmetrical 
points within the Brillouin zoneJ^ The Mermin- Wagner 
theorem is not valid for two-dimensional magnets with 
a dipolar coupling of spins with continuous degeneracy, 
and a true long range order can exist even for a purely 
two-dimensional case at finite temperatures*^ - — 

The models of dipolar magnets were discussed 
originally in regard to real crystalline spin systems. 
During the last decade the most impressive achievements 
in magnetism were related to fabrication, investigation, 
and application of artificial magnetic materials, see 
Ref. 120 for a recent review. The nanotechnologies today 
have progressed to the state where the manufacture 
of nanosize, periodic magnetic superlattices of different 
types is feasible. Among them two-dimensional lattices of 
sub-micron magnetic particles (so-called magnetic dots) 
attract much attention. These magnetic dots, of different 
forms and of a submicron size, are made of soft magnetic 
materials such as Fe, Ni, Co and permalloy^ - — or highly 
anisotropic materials like dysprosium^ or FePti21 



In the dot array lattice dots are separated from 
each other so that direct exchange interaction between 
the dots is negligible. Thus the dipolar interaction 
is the sole source of coupling between dots and the 
configuration of dot magnetic moments is dictated by 
the dipolar interactions of the dots and by the external 
field. Owing to the absence of exchange, magnetic dot 
arrays constitute promising material for high-density 
magnetic storage media. For this purposes, the dense 
arrays of small enough magnetic dots with the magnetic 
moments perpendicular to the array plane are optimal, 
see Refs. 28-30, Currently, ordered arrays of magnetic 
submicron elements have been discussed as materials 
for so-called magnonics, i.e., a new field in the applied 
physics of magnetism in which magnon modes with a 
discrete spectrum present for magnetic nanoelements are 
used in devices for processing microwave signals. 31 

For small enough dots with size of the order of 100 
nm the magnetization inside of a dot is almost uniform, 
producing the total magnetic moment mo 3> ^b, where 
fi b is the Bohr magneton, the typical value for an atomic 
magnetic moment. For rather small magnetic dots of 
volume 10 4 -10 6 nm the value of mo exceeds 10 /Ub, and 
for dense arrays the characteristic energy is higher than 
the energy of thermal motion at room temperaturej 1 1 i 12 
The individual dots in an array do not touch each 
other, and their interaction is only determined by the 
dipole interaction of the magnetic moments. Therefore, 
magnetic dot arrays represent a new kind of magnetic 
material with purely two-dimensional lattice structure 
and high enough pure dipolar coupling between magnetic 
moments. Such systems represent dipolar magnets and 
fill their theoretical investigation with a new physical 
content. Thus, magnetic dot arrays are interesting as 
radically new objects for the fundamental physics of 
magnetism. 

We will discuss only the situation where the magnetic 
moment of an individual magnetic dot is perpendicular 
to the array plane (cn/-plane) , m = ±moez, and the 
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system can be described on the basis of the Ising 
model. This situation is most promising for perpendicular 
magnetic recording systems. The energy of dipolar 
interaction of Ising moments perpendicular to the 
system's plane is minimal for antiparallel orientation 
of magnetic moments, that, of course, can not be 
fulfilled for any pair of particles in the array. Within 
the nearest-neighbors approximation, such interactions 
lead to antiferromagnetic (AFM) structures, e.g., simple 
chessboard AFM ordering is known for two-dimensional 
Ising square lattice with dipolar interaction^ To explain 
the validity of this result for real dipolar interaction it 
is enough to mention that accounting for next nearest 
neighbors gives a correction to the energy of the order of 
30 % only. 

Thus, the magnetic structure for square lattice dot 
arrays can be easily understood. However, the close- 
packed triangular lattices of the magnetic dots are 
also frequently used in experiments. In particular, these 
lattices of cylindrical particles considerably extended in 
the direction normal to the array plane are naturally 
obtained when the array is prepared by controlled self- 
organizatiom 3 - 2 - Again, the properties of these systems can 
be described in the two-dimensional Ising model with 
AFM interactions. However the triangular lattice with 
AFM interaction of the moments is a typical example of 
frustrated antiferromagnets, see for review ^ 

For frustrated magnets spins interact through 
competing exchange interactions that cannot be 
simultaneously satisfied, giving rise to a large degeneracy 
of the ground state of the system. For nearest- 
neighbor Ising triangular lattice with AFM interaction, 
the thermodynamic properties are quite unusual^ It 
is enough to mention that in this model magnetic 
ordering is absent at any finite temperature T ^ 
0; the ordering appears as a result of accounting 
for the next-nearest-neighbor interactions only^ This 
counterintuitive feature can be explained within the 
concept of creation of linear topological defects^. Thus, 
in contrast to bipartite square lattice, nearest-neighbors 
approximation did not provide even adequate zero 
approximation to the problem of the ground state of a 
triangular lattice. 

In the present work the ground state of a triangular 
lattice of mesoscopic magnetic dots, each having a 
strong easy axis for magnetization perpendicular to the 
array plane, and in an external magnetic field also 
perpendicular to the plane of the dot lattice, will be 
considered. Both unbounded planar lattice and various 
bordered semi-infinite or finite elements of the triangular 
lattice are investigated. A cascade of phases with different 
patterns of dot magnetization has been found; these 
constitute the sequence of ground states as a function 
of the external magnetic field. In contrast to a square 
lattice, the transition between these states is governed by 
a novel mechanism involving creation of linear topological 
defects with non-zero magnetization^ 



II. MODEL DESCRIPTION. 

Consider the set of Ising magnetic moments m„ = 
<7TOo e zj f i 1) each parallel to the z-axis, and placed 
in a sites of a triangular lattice n, 

n = ake x + y (e x + VSe.yJ . (1) 

where a is a lattice constant, fc, / are integers, and e x and 
e y are unit vectors parallel to x and y axis, respectively. 
The magnetic moments are interacting through the 
magnetic dipole interaction, and an external magnetic 
field H = He z is applied perpendicularly to the array's 
plane. The Hamiltonian of this system of magnetic 
moments can be written as 

^^oE-^p-offE^, (2) 

where the first term describes dipolar interaction, with 
the summation performed over all of the pairs of the 
lattice sites. Below for the sake of simplicity we will 
present the energy (per one magnetic particle) in the 
units of ra§/a 3 and we will use the dimensionless 
magnetic field, h — H/H*, where characteristic value 
H* = mo /a 3 . The present model is clearly not restricted 
to a dot lattice of the type explicitly described above, but 
also applies directly to any triangular lattice of identical 
dipoles that are restricted to the two directions of normal 
orientation^ It is interesting that the model formulated 
in this paper can be used to describe a system of a vortex 
state magnetic dots, accounting for the interaction of a 
magnetic moment of vortex cores ^ 

As has been mentioned above, the triangular 
lattice is known as a typical frustrated lattice for 
antiferromagnetic ordering, and the frustration is present 
even for simplest nearest-neighbor interaction. Of course, 
not only nearest neighbors are important for magnetic 
dipole interaction. The dipole-dipole interaction is long 
ranged, that frequently leads to quite complicated 
structures with many sublattices and with high level 
of degeneracy. The presence of these two sources of 
degeneracy makes the problem less definite, and it is 
not obvious a priori what structure will constitute the 
ground state for such a lattice. In this situation it 
is natural to start with the numerical analysis of the 
problem. 



III. GROUND STATES FOR INFINITE 
SYSTEM: NUMERICAL ANALYSIS 

We perform Monte-Carlo analysis (simulated 
annealing, see for details Appendix) of the magnetic 
configurations with minimal energy at zero magnetic 
field as well as for for different values of the magnetic 
field. 



Simplest ground states of the system: zero field 
and saturation. 



It is clear that for high enough magnetic field all 
magnetic moments will be parallel to the field, giving 
the saturated magnetic structure, which can be referred 
to as a ferromagnetic structure. The ground state of the 
array in the absence of magnetic field is much less trivial. 

It is known that for triangular Ising lattice in the 
several-neighbor approximation the simple AFM order 
with two sublattices can be implemented^ As well, we 
found the same configuration for the long-range dipolar 
interaction in the absence of the field and for a small 
enough magnetic field. For these states, the magnetic 
elementary cell is rectangular having lower symmetry 
than for the underlying triangular lattice, and this state 
possesses much higher discrete degeneracy than the 
simple chessboard structure for a square lattice discussed 
before. Several such AFM states can occur in the 
system, which are different but fully equivalent in their 
energies. Fig. Q] presents three of these states, while the 
other three states are obtained from them by changing 
the magnetic moment sign <r n at all of the particles. 
Then, we will briefly describe the magnetic state of 
a given particle, specifying the number of neighbors 
with the favorable and unfavorable orientations. The 
antiferromagnetic states present on the figure Q] can be 
called the 4-2 type state. Of all of the states of the system 
with the simple periodic distribution of the moments, 
this 4-2 type state corresponds to the minimal energy. 
It is not optimal because only four neighbors of the six 
nearest neighbors of each moment have the sign opposite 
to the sign of this moment, while the remaining two 
neighbors are parallel, and this is unfavorable to the 
AFM interaction. The necessity of deviation from the 6-0 
optimal structure is a typical manifestation of frustration 
in the system. 
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Phc. 2: Mean value of the magnetization (mo) (in units of 
mo, per one dot) of the array as a function of magnetic field 
(in units of H, = Mo /a 3 ) found by Monte-Carlo simulations. 
Detailed data at low fields are present at the insert. 



magnetic moment is present. Monte- Carlo analysis shows 
that for some small but finite values of the magnetic field, 
at least up to h — 0.7, the mean value of the magnetic 
moment (to) equals zero indicating to the simple AFM 
structure. For higher fields, numerous more complex 
structures with < (m) < too occur in the intermediate 
region between AFM state and saturated state. The mean 
value of the magnetic moment (per one particle) (m) 
corresponding to these configurations, is present on the 

Fig.m 

Note the specific regions of this dependence, present at 
different field intervals; first, the region with small values 
of (to) < 0.2too having rather non-regular dependence 
of (to) on h; second, the regions with constant values 
of (m) independent on the magnetic field (shelves); and 
third, the saturation region. The characteristic magnetic 
structures found in these regions are depicted at the 
Fig. SI 

Monte-Carlo data are not too clear in the region of 
small fields such as 0.7-0.9, and the magnetic structures 



are far from the simple AFM structures, see Figs. 3(a) 



[3(d)] and we will discuss the magnetic structure within 
this region below. 



Phc. 1: Uniform states of the 4-2 type giving the energy 
minimum at small magnetic field. Here and below at all figures 
the open and closed circles denote the particles with the 
upward and downward moments, respectively. 



B. Monte-Carlo analysis for intermediate field 
values. 



As was mentioned above, for a zero magnetic field 
the simple AFM structure with zero mean value of the 



1. Shelves and saturation. 

Within the shelf regions, almost all initial Monte-Carlo 
configurations lead to the same magnetic structures, 
which correspond to the formation of triangular 
superlattices for the minority dots (antiparallel to the 
magnetic field) with different lattice spacings. As an 
example, note the ideal triangular superlattices with 
(to) = too/ 3 and with the period a s \/a = -\/3, see fig. |3(f)| 
present at the values 2.4 < h < 6.4. For h igher fields, the 
superlattices with a s \/a = 2 (Fig. |3(h)[ ), a s \/a = \fl ■, 
(Fig. 3(i) ) and a s \/a = 3 in Fig. 3(j)[ ) correspond to such 
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(c) h = 1.5 
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(b) h = 1.1 



o«o om 
o»«oo««o 

•o o«o 

OMOO«» 

o«o o«oo 
•oottooi 

ooto omo 
' ' oot#oo 



(d) h = 2.1 
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(f) h = 4.0 



oooooooooo 
• • • • • 

oooooooooo 
• :• • • • 

oooooooooo 
•o»o»o»o»o 
oooooooooo 
•o»o»o»o»o 

r . OOOO 

•o«o«o«o»o 



(h) h = 7.5 
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for the shelf regions. For example, at the values 0.9 < 
h < 1.5, where the finite (but small) magnetic moment 
(mo) is formed, the translational symmetry for the set of 
flipped dots cannot be attributed to simple superlattice 
structure, see Fig. 3(a) But in this figure a novel 
element, the additional zigzag line of the sites, oriented 
parallel to one of the translation vector of the lattice, 
is clearly seen. The resulting magnetic structure can be 
interpreted as an antiferromagnetic domain structure in 
the system, with the zigzag line as a domain wall. Such 
a scenario, magnetization through the creation of a set 
of topological linear defects, was described for a two- 
sublattice antiferromagnetic state with an interaction of a 
few neighboring moments. 37 For the region of small fields, 
the increasing the magnetic field leads to an increase of 
the density of the topological linear defects, see Figs. |3(b)| 



and Fig. 3(c) 



The common "topological" scenarios are present for 
other transition regions; both below and above the 
shelf regions with the magnetic structure of a form 
of ideal triangular superlattices of minority dots. For 
example, the structure present at 2.0 < h < 2.4 
can be described as a "compression" of the domains 
of the superlattice of period ay/3 by the lines of dots 
with down magnetic moments, see Figs. |3(d)| and |3(e)| 
whereas the state at the opposite end of this shelf can 
be seen as a "rarefication" of the a\/3 superlattice, see. 
Fig. 3(g) The transition structures corresponding to the 



(j) h = 10.1 



"higher" shelves have Cq symmetry, higher than for low 
field structures. Note the essential difference of such 
topological behavior with what is known for a square 
lattice, where the competition of square and triangular 
(distorted) superlattices of minority of magnetic dots is 
responsible for magnetization processes. 



C. Small fields; exhaustive search of the states. 



Phc. 3: Ground states for characteristic values of magnetic 
field found by Monte-Carlo simulations. 



shelves. 

For the high magnetic fields near the saturation region 
the magnetization process is going through the same 
scenario as for a square lattice, namely, by a flip of a small 
amount of magnetic moments and creating a superlattice 
of flipped dots of small density. 



2. Transition region, topological mechanism. 

In the region of low magnetic fields, as well as in the 
regions of magnetic field where the transitions between 
the superlattices occur, resettability of Monte-Carlo 
result is lowering, and the results becomes unreliable. 
The observed magnetic structures in these transition 
regions are characterized by much lower symmetry than 



As has been found by direct Monte-Carlo simulations, 
in the field region 0.9 < h < 1.5 the ground state is 
realized by a system of parallel AFM stripe domains of 
the wight depending of the field, see Figs. 3(a) - 3(c) The 



minimal field for the start of this process corresponds to 
a low density of such defects, and to find the critical field 
one needs to consider larger and larger system. Namely, 
to present a stripe of weight n we need a system of size at 
least (2n+l) x (2n+l). Below in the Sec.|IV]we will find 
the starting field for the creation of the set of topological 
defects by an analytical calculation. 

For refinement the Monte-Carlo data and for 
clarification of the magnetic states at the fields of 
interest, < h < 3, we perform the direct exhaustive 
search based on the picture of stripe AFM domain 
structures for rhombus-shaped space regions with various 
(not necessary equidistant) geometries of domain lines, 
up to the size 60 x 60. It appears that for all fields the 
only equidistant structures corresponds to the minimal 
configurations, with linear system of stripes at h < 1.5 
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<m> 




Phc. 4: Magnetization function at low fields found by 
exhaustive search of the states of rhombus-shaped samples 
(full line) and the Monte-Carlo data (symbols). 



or triangular superlattice at h > 2.0. Then the only 
equidistant structures with the size up to 300 x 300 
were examined. The magnetization curve based on these 
calculations is represented in the Fig. 2] and compared 
with the Monte-Carlo data. 



IV. MAGNETIC GROUND STATES: 
DESCRIPTION. 



ANALYTIC 



As has been found in the previous section, for small 
values of the magnetic field, the magnetic moment 
equals zero corresponding to the simple two-sublattice 
AFM structure, whereas for high magnetic field the 
saturated magnetic structure, or ferromagnetic structure, 
is present. It is known that for a non- frustrated square 
lattice of Ising magnetic moments the destruction of both 
types of the magnetic order is started through creation of 
the point defect in the state, the single magnetic dot with 
the magnetic moment reversed with respect to the regular 
structure of a given stated First, in this section we will 
check the validity of such scenarios for ferromagnetic and 
AFM states of triangular lattice array of magnetic dots. 
Then the theoretical description of the novel topological 
mechanisms will be obtained. 



A. Point defect scenarios. 

In order to determine the values of the magnetic fields 
which correspond to such "point defect instability" we 
have calculated the change in dipolar interaction energy 
that occurs when the magnetic moment of a single 
dot is reversed with respect to the ferromagnetic and 
two-sublattice AFM structures. This energy change is 
determined by the energy per dot in the initial states^ 
which can be expressed by simple lattice sums calculated 
with high precision. These sums here and below were 
calculated with standard program package Mathematica. 



1. Saturation. 

The point defect scenario describes well the instability 
of the saturated state. It is easy to see that the change of 
the energy of the saturated state with the flip of a single 
magnetic moment can be presented as 



W x = 2m {H - H s: 



(3) 



where 



ffsat = m T^TE = h-t^T ' h ^ ~ U-034176 . (4) 

n^O |ni1 a 

The quantity i/ sa t determines the saturation field for 
a triangular lattice of Ising magnetic moments. If the 
magnetic field H < H sat , the value of W\ is negative 
and the flipping of a dot becomes favorable. The value 
h s& t is higher than for a square lattice, /i^" 3,10 = 9.034, 
which just reflects the higher density of the triangular 
lattice comparing with the square lattice. Note here that 
the accounting for only nearest-neighbor interaction gives 
a much lower value = 6 clearly demonstrating the 
importance of the long range character of magnetic dipole 
interaction for the description of the ferromagnetic state. 

The energy of the simultaneous flip of a pair of 
magnetic moments placed at the distance \n\ can be 
easily written as 



W 2 = 2W 1 +4m 2 /\n\ s 



(5) 



that can be interpreted as a repulsion of the flipped dots. 

Evidently, for H < H sat some magnetic moments tend 
to reverse, and the creation of a finite density of reversed 
magnetic moments becomes favorable. The flipped dots 
should be dispersed as far from each other as possible, 
in order to minimize the energy of their repulsion. Thus 
one would expect these flipped dots to be organized in 
a triangular superlattice, with the lattice spacing a s \. 
We can regard the resulting system as the superposition, 
on the original triangular lattice with magnetic moment 
+mo, of the triangular superlattice of "double dots" 
of moment magnetic moment — 2mo- From this simple 
picture, it is clear that the mean value of the magnetic 
moment per one dot (to) can be expressed through a s \ as 
following, (to) = mo — 2mo(a/a s i) 2 . 

This picture also allows us to present the energy of 
dipolar interaction of the array with a small density 
superlattice through the known value of h sa t , see Ref.njl. 
First, note that the flipped (double) dots experience 
the field i/sat from the rest of array. Consequently 
the interaction energy of the triangular lattice and 
superlattice is 2moei? S at, £ = (mo — (to)) /mo, per one 
dot of the full lattice. Then the contribution to the 
energy of each double dot due to its interaction with 
all other double dots equals to (l/2)/i sat (2TOo) 2 /(a s i) 3 , 
per one double dot. This contributes AWfm (a/asi) 5 ^ 2 
to the mean dipolar energy per one dot of the overall 
lattice. Combining all these contributions and adding the 



Zeeman energy Wh, Wh = —{m)H, we can obtain the 
approximate formula for the energy of the superlattice 
state per one dot in the form 



,2 r 



W = 



''sat 



-h + e(h- hsat) + 2h sat (|) 



5/2* 



Minimizing the energy © over (m) (in fact, over e) yields 

'2(h sat -h)\ 2/3 



= 1-2 



ra 



5/i s 



(6) 
slds 

(7) 



This dependence describes very well the numerical 
data for the dependence (m) on H near the saturation, 
which can be fitted by the dependence common to (J7J), 
with the power 0.656 instead of 2/3 ~ 0.667, and with the 
coefficient 1.087 instead of 2(2/5) 2 / 3 ~ 1.086. Thus the 
point defect scenario based on the reversal of a single 
magnetic moment describes well the instability of the 
saturated state for both triangular and square lattices. 



2. Low fields. 

In principle, the common calculations can be 
performed for AFM state, as well as for any state with the 
superlattice of flipped dots of the same symmetry as for 
underlying dot lattice. This approach describes well the 
instability point for the AFM state for a square lattice of 
magnetic dots^i Let apply it to our triangular dot lattice. 
Reversing the magnetic moment of one dot in the AFM 
state becomes favorable at H > -Hafm = 2-Eafm/to-o, 
where 



FM 



mo \ ^ 

~ a 3 ^ 

k,l^0 



mo 



[3 (2fc + l) 2 + (2Z + l) 2 

E 



3/2 



[3fc 2 + P] 3/2 



1.8377 



to 



(8) 



This value is much higher than the instability field h ~ 
0.7 found numerically in the previous section. Thus the 
reversal of a single magnetic moment cannot describe the 
instability of the AFM state for triangular lattice. 



B. Instability of AFM state through creation of 
topological line. 

As we found by Monte-Carlo analysis, the AFM 
state looses its stability as the field increases because 
of the creation of topological defect lines. These lines 
can be also called domain walls, because within the 
description of AFM structure in terms of the standard 
antiferromagnetic order parameter (antiferromagnetism 
vector L = m.! — m 2 , where ni! and m 2 are the magnetic 
moments for different sublattices) the values of L have 
different signs on both sides of this wall. This defect line 



corresponds to an additional zigzag line of the particles 
with magnetic moments, which are oriented as their 
neighbors, and should be normal to one of the elementary 
translation vectors, see Fig.0 
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Phc. 5: (color online), (a) - (c), the structures for topological 
defect; (d) the ideal AFM structure; the characteristic lines 
used for calculation of the defect energy are shown, see details 
in the text. As usually, the open and closed circles denote the 
particles with the upward and downward magnetic moments, 
respectively, but the dots with upward magnetic moments 
within the defect line are mentioned by grey. 

An importance of defect lines for thermodynamics is 
a well-known property of two-dimensional systems with 
discrete symmetry breaking. Because of the creation 
of a finite density of such lines, the long range order 
is destroyed at finite temperature, determined by the 
energy of the defect, see, for example, the articled But 
for frustrated AFM states the behavior can be very 
unusual. In particular, the magnetic order for AFM Ising 
system with the nearest-neighbors interaction is absent 
at any finite temperature T > . 33 i 34 Generally, this 
behavior can be explained using the defect line picture of 
the phase transition, with the vanishing of energy (more 
exactly, the free energy) of a certain linear topological 
defect. For triangular lattice Ising model with nearest- 
neighbor interaction it could be linear twin boundary 
defects, which separate the different states presented 
in Fig. [TJ 35 Examples of these states appear for finite 
systems, see below Sec. I VII and figures wherein. The 
formation of these defects does not violate the 4-2 
relation. Thus their energy is zero in the nearest-neighbor 
interaction approximation and should be positive but 
small for dipolar coupling. 

For our system of mesoscopic magnetic particles with 
a dipole interaction, the effect of a magnetic field, instead 
of thermal effects, should be significant. The topological 
defect line with non-zero magnetization was recently 
found. 37 This defect coincides with that observed in our 
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numerical simulations, compare Fig. [S] and Figs. 3 (a 



|3(d)| above. The particles directly entering into the 
defect line have an unfavorable configuration of the 3-3 
type, and their presence results in energy loss; however, 
they are adjacent to the particles (with the downward 
magnetic moment) having a 5-1 configuration more 
advantageous than the standard 4-2 configuration. The 
numbers of these anomalous states coincide with each 
other. Thus the creation of this wall does not result 
in energy loss in the system in the nearest-neighbor 
approximation, and the loss should be small for particles 
with dipolar coupling. In contrast to twin boundary 
defects considered before^ this domain wall has a 
nonzero magnetic moment. As a physical consequences 
of this property note that the behavior of the system in 
a magnetic field is dictated by the defects with nonzero 
magnetic moment. In fact, for such a defect the additional 
energy gain m$H per defect particle appears in the 
magnetic field H . Then the defect energy decreases as 
the field increases and becomes zero at H = i/DW = 
Euw/ m o- Then for H > -/few the finite density of such 
defects will be present in the ground state. This is exactly 
the scenario observed in our numerical simulations at 
magnetic field at the range 0.7 - 1.5. 

In order to find the critical value of the field ffow, 
let calculate the energy of the domain wall £dw- It is 
convenient to use the symmetry of the state with the 
defect line and to divide the full lattice into lines of dots 
parallel to the defect line, as it is shown in the Fig. [5] 
It is evident that all dots within one such line have the 
same energy, and the energy of the domain wall (per one 
dot in the defect line) can be present as a sum over these 
lines as following 



few — ~7J ^ m QV n H n 



(9) 



where the integer n describes the distance a n of the given 
line from the defect line, a n — an/2, a n — ±1 gives 
the sign of the moment for the n-th line, and H n is the 
magnetic field created on the dot in the n-th line by other 
dots in the system. To find the field H n it is convenient 
to group all other dots to pares of lines, equidistant from 
the n-th line, as it is shown for n = 0, 1, 2 at Fig. (a), 
(b) and (c), respectively. Let us enumerate these pairs by 
an integer k so that distance between n-th line and one 
component of the fc-th pair equals to afc/2, the pairs with 
k = 1,2,3 are present at the Fig. [5] Then the energy of 
the magnetic state with a domain wall can be presented 
by a double sum, over n > and k > 0. 

It is easy to see that for any finite n the only pairs 
with limited k < n contribute to the energy of the 
state with the domain wall. For example, for the lines 
directly entering the defect line [n = 0, see Fig. [5] (a)] 
the contributions of two lines composed any pair cancel 
each other. For this line, the non-zero contribution to the 
energy is given by the dots from the same line, denote 
this contribution as £o. Then, for the line with n = 1, 
only one pair gives non-zero contribution, see Fig. [5] (b) , 



and the energy can be written as £o — 2s\. Similarly, for 
n = 2 the energy is £o — 2ei + 2e 2 , see Fig. [5] (c), and 
so on. Finally, the energy of the state with domain wall 
is presented through Eq and the particular finite sums of 
the positive quantities e n , 



£2n+l — ^ 



e 2n - 



£j[(„ + 1/2)3 + 3(fc-l/2)*]3/a 
4 2 



^ (n 2 + 3k 2 ) 3 / 2 n 3 



(10) 



The energy of the domain wall equals the difference 
of the energy of the state with the domain wall and the 
ground state energy. To find the ground state energy, it 
is convenient to use the same presentation by the parallel 
lines, see Fig. [3] (d), and to present it by the same sums 
e n . It is clear that the energy per one dot in any line 
in the ground state is proportional to an infinite sum of 
the form £qs = £o + 2 X)^Li(~l)" £ n- Then the domain 
wall energy can be found by term-by-term summation 
of corresponding contributions of the form [(eq — £gs) + 
(e - 2ei - e G s) + •••] = few = 2ei - 4e 2 + 6e 3 + .... The 
corresponding infinite series few = — 2 Y^n=i ( — l)™ ne n 
are sign-alternating and converges quite well. Finally, 
domain wall energy per one dot -few can be presented 
as follows 



E 



DW 



MO^DW , ^DW — few 



m 



h 



DW 



0.70858944 . 

(11) 

Here we also present the characteristic value of the 
magnetic field, and -/few = £>dw/t*o determining the 
border of stability of the simple AFM state; for H > 
ifew AFM state becomes unstable against creation of 
domain walls. Note the calculated value (fTTj) is in good 
agreement with that found by numerical simulations, but 
it is much lower than the field of point defect instability 
for AFM state, H AFM = 1.8377m /a 3 . 



C. Plateau description 

Monte-Carlo simulations result in some peculiarities 
(See Fig. [6|) in the dependence of the magnetization on 
applied magnetic field in the form of plateaus, where the 
value of the function does not change over a wide range 
of the argument. 

These peculiarities have a simple explanation. The 
magnetization of the array increases at small external 
field due to the formation of parallel topological 
defects in the form of domain walls. At some critical 
concentration of such walls the resulting state in 
nothing but the superlattice of flipped dots which 
has the triangular structure that coincides with the 
array symmetry; See Fig. [3^f,h,i,j). Such a superlattice 
transforms into the structure similar to itself but 
with the other step (lattice constant) as the applied 
field increases. Since the lattice constant has discrete 
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Phc. 6: (color online) The magnetization function dependence 
on external magnetic field. Numerical data are present by 
solid line, analytical calculations of the instability points 
specified by the symbols, see Eq. I13I 



values a sl /a = y/3,2, V7,3,2\/3, Vl3,4,\/l9, V21,5..., 
then such superstructure has good stability against the 
alteration of the external field and magnetization can be 
changed only stepwise. One can see that such a structure 
consists of two inversely magnetized ferromagnetic states 
with the lattice constants a and a s \. As stated above the 
magnetization of a such state is 



(to) = m - 2m a (a/a s \Y 



(12) 



The value of the field of the stability loss of such 
superstructure relative to the transition to other lattice 
constant can be easily calculated on the same principle 
as the field of the ferromagnetic state stability that was 
done above. 



mo 



(13) 



where the multiplier 2 in the numerator responds to the 
change of magnetization of the dot in comparison with 
the ferromagnetic state and h sat — 11.034176 is the field 
of the transition to the saturated state. The values of 
these fields depending on the superlattice magnetization 
are represented by symbols (stars) on the Fig. |B1 

As was found by numerical simulations, the destruction 
of such superstructures can go through the creation 
the lines of topological defects, see Fig. [3jg) . This 
mechanism assumes that unit cells of the superlattice 
repel themselves with an increase of the field and the 
line of the dots magnetized inversely to the dots of the 
superlattice passes between them. The situation recurs 
until a new superlattice formes with a larger lattice 
constant. Though this effect occurs in the narrow range 
of the field it leads to instability of the superlattice at 



smaller value of the field than is predicted in equation 

<tn- 

Formation of the triangular superlattice is impossible 
at small values of the applied magnetic field h < 2.5 
but some rectangular superlattices with C2-symmetry are 
possible, see Fig. GH[c), which also are the reason for the 
plateaus appearance (but not so well pronounced) at the 
magnetization function, see Fig. @] 



V. SEMI-INFINITE ARRAYS: THE ROLE OF 
THE BOUNDARIES. 

In the previous theoretical consideration we analyze 
the idealized model of the infinite array. Of course, real 
superstructures are large, but finite systems. The border 
elements (edge surface of an array) are expected to play 
a considerable role in the formation of the properties 
of the ground state. The existence of any translational 
invariance significantly simplifies the problem. For finite 
system, there is no translational symmetry. The finiteness 
of the system manifests itself in two different ways; first, 
through the difference in the coordination numbers for 
the dot at the border and in the bulk; and second, by the 
direct influence of the system size. In this sections, we 
will analyze the "border problem" only, assuming that 
the array is semi-infinite, e.g., it is bordered by one 
border line. As well we will discuss the "edge problem", 
investigating the properties of an edge element, a single 
dot located at the crossing of two borders. The systems 
of finite size will be considered in the next section, 
mainly numerically. We primarily will discuss only the 
simplest geometry of the system, supposing that the 
array's borders are parallel to simplest translation vectors 
of the underlying triangular lattice. 



A. Analytical description 

Consider the simple two-sublattice AFM state typical 
for an array at small enough magnetic fields. For an 
infinite array, any magnetic dot is influenced by the 
magnetic field, generated by other dots of the array. 
This field is parallel to the dot's magnetic moment; 
and the strength of the field is given by the expression 
^■afm = 1.8377. 

Let consider some different kinds of bordered arrays. It 
is clear that the same quantity, the dipole magnetic field 
at a dot located near the border, determines the stability 
of the state. In order to calculate the field on a given 
dot, it is convenient to take the coordinate system with 
the origin at this dot, see Fig. [71 For all cases of interest, 
namely, for a magnetic dot located at the edge of the 
array or at one of its borders of different orientation, a 
certain common property is easily seen. The magnetic 
field acting at the dot can be expressed in terms of 
auxiliary sums, from the single sums over dots located 
at the ray beginning at the coordinate origin (half the 
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coordinate axis) and a few double sums over dots located 
in one of the array sectors that is shown in Fig. [7J Then 
for concrete estimates we will chose the dot at the origin 
with the "down" magnetic moment, so a negative value of 
the field corresponds to stability of the configuration and 
a positive value of the field corresponds to the instability 
of the state. 




Phc. 7: (color online) The definition of auxiliary sums used 
for analysis of semiinfinite arrays in AFM state, see the text. 

The single sums are expressed in terms of the Riemann 
C(3)- function, 



oo ^ 

?i = V — = C(3) = 1-20206, 

Z ✓ pO 

71=1 
OO 

^ = E 



^— '— = -3C(3)/4 = -0.90154, (14) 

<K1 <J 



1 C(3) 
^ h(3n*f 2 3^3 



0.23134, 



whereas the double sums can be easily determined 
numerically, for example, 



asym 

T 60 — 



oo oo 



EE 



(-1)" 



n=i k=i K + kn + k 



21 3/ 2 



-0.22164 (15) 



<C = ~ [hAFM + ^Z Vm + + 2Ci] - 0.12545. 

(16) 

Using a simple geometrical consideration, all the double 
sums of interest can expressed through the value of 
h-AFM = 1.8377 and the partial sum Cgg^™ found 
above. The knowledge of this particular sums gives the 
possibility to determine the stability of any particular dot 
of interest. 

Let us start with the infinite borders parallel to the 
elementary translation vector. There are two kinds of 
such a border, with all magnetic moments parallel and 
with the alternating magnetic moments, see Fig. [7] The 



magnetic field on the dot in these borders, and hf±, 
respectively, can be written as 



fi - h AFM /2 = 0.2832, 



C-i - ^afm/2 - -1.8204. (17) 



Note first that is positive, and the infinite borders 
with parallel magnetic moments is unstable in the 
absence of the magnetic field, such a fragment can appear 
at some finite value of the magnetic field only. Both 
the aforementioned properties are in agreement with 
the numerical analysis of finite arrays, see the next 
section. In particular, the border with parallel magnetic 
moments never appears for finite arrays; instead, the 
complicated multi-domain AFM structure is present. 
It is worth noting an essential difference in behavior 
between triangular and square lattices. In the latter case, 
the simple AFM state with ^-symmetry (chessboard 
AFM) is the ground state for finite systems with any 
shapes of the array including the systems with acute 
angles or for system of circular shape. 12 . In contrast, 
the "perfect" AFM ordering with (^-symmetry for the 
triangular lattice is possible only for rhombic shape of the 
sample, borders of which are parallel to the fundamental 
translation vectors of the lattice. Lines which consist of 
unidirectional dots would have sharp bends in all other 
cases, See Section fVTl 

The next point of interest is the behavior of the 
magnetic dot with the downward magnetic moment at 
the edge of the array in the AFM state. For a square 
lattice, the magnetic moment of such dots becomes 
unstable at the magnetic field that is much weaker than 
the instability field for the infinite system, h s J^ re = 
1.563. Further, the magnetic moment of this dot is 
reversed at this value of the field. For a square lattice, 
such reversal of the edge dot is the beginning of the 
destruction of the AFM state in the finite array. For 
the triangular lattice, the edges of finite array are also 
a "weak points" for a single-dot instability. The values of 
the instability field for a dot at the vertex of edges with 
the angles 60° and 120°, can be expressed through the 
sums <7gQ m and 0120, respectively. The state of the edge 
dot is then stable and does not change until the external 
field increases to the value. The field h^o of the reversal 
of the single dot at the vertex of 60-degree edge can be 
present as 



/ sym 

w 



'(ill 



■2£_i) = 1.67764, 



(18) 



whereas for 120-degree edge the instability field /1120 has 
the lower value, 

h 120 = -(2<a Vm + £1 + 2£_i) = 1.0443 (19) 

Thus we found that the more sharp edge appears to 
be more stable against the action of the magnetic field 
(note that for square lattice the situation is opposite^). 
This result seems to be a contra-intuitive, but it is 
the reflection of the frustration for AFM state at the 
triangular lattice. 
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The above field for single-dot edge and border 
instability are weaker than the instability field for the 
border of infinite system hfi = 1.82039 but their values 
are higher than the value hnw = 0.7085, where the 
creation of the domain wall becomes favorable. Thus, 
the scenario with the creation of such domain walls is 
preferable in a thermodynamic limit. On the other hand, 
the creation of a topological domain wall require the 
perturbation with a reversal of an essential number of 
dots (formally, proportional to the system area). Thus 
such a process needs the overcoming of a very high 
(formally, infinite) potential barrier, and it is hard to 
realize it for a real experiment. 

To find an alternative scenario, let us consider the 
domain wall which pass nearby the 120-degree corner. 
Thereto we fix the corner site to the origin of coordinates 
and build y-axis on the bisector of the angle, see insert 
on the Fig. [5] Then the domain wall which lays at the 
distance i/dw results in the magnetization flip-over, in 
comparison with the ground state, on all sites left from 
the wall. Single dot on the vertex angle is equivalent 
to the domain wall with the coordinate yow = 0. We 
calculated the necessary field of such a domain wall 
appearance for different i/dw, see Fig. [5J Calculations 
were carried out by the direct summation of the array 
energy, which consists of 400 lines in the y direction. 

As one can see from the the data on Fig. [SJ the 
necessary field is higher for the domain wall near the 
corner of the array than for the deep-laid wall. Also, 
we can see from the figure that the field of the domain 
wall generation depends on the relative directions of 
the applied magnetic field and the magnetization of the 
vertex site (taken, as before, negative on the figure). 
Thus only domain walls with integer values of i/dw are 
of interest for the problem. With an increase of yow 
the field that is required for the wall formation decreases 
converging exponentially to the field of creation of infinite 
wall, that is similar to the presence of the potential 
barrier on the surfaces of solid. Therefore one can assert 
that the smaller the difference between applied field and 
how = 0.7085 the farther from the corner domain wall 
appears. 

It is necessary to take into account for large arrays 
that domain wall at the sufficient distance \y\ from the 
site produces the magnetic field at this site proportional 
to 



dx 



(x 2 + y 2 ) 



2^/2 



y 



(20) 



notably converges to zero when \y\ — > oo. Therefore, 
the formation of the another deep lying domain wall 
is more energetically favorable than the wall near the 
corner for the small concentration of topological defects. 
And consequently the magnetization of the large arrays 
increases firstly in the volume, where equilibrium density 
of parallel topological defects would be observed; and 
only in the case of external field considerably exceeding 
h-DW = 0.7085 domain walls go up to the corner. This 
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Phc. 8: Dependence of the field of line formation on the 
distance to the corner of the array. 



case is opposite the magnetization process for a square 
lattice which is initialized on the borders. If we go out of 
the thermodynamical limit and look at the slow dynamics 
of the system then we could expect that for the strong 
external magnetic field (close to hvzo) the topological 
defects initially appears close to the corner and then 
drifts to the volume due to the aforesaid reduction of 
the potential barrier. 



VI. GROUND STATE FOR FINITE ARRAYS 

An analysis of finite arrays was accomplished 
numerically, with usage of the same Monte- Carlo 
procedure as described above for infinite systems. The 
different shapes, either fully consistent with the geometry 
of triangular lattice (bordered by dense lines of magnetic 
dots only, like rhombus, triangle, hexagon) or less 
consistent with underlying lattice (rectangle, circle) 
shows a big variety of the magnetic structures. 



A. Zero magnetic field 

As has been found in the previous section, the 
simple two-sublattice AFM survives near the border 
of semi-infinite system, with one essential exception: 
the border of the type of a dense line with parallel 
magnetic moments is unstable at zero magnetic field. 
This instability is crucial for an understanding of the 
formation of the ground state for finite samples bordered 
by dense lines only (rhombus, triangle and hexagon). 



1. rhombus 

The rhombus is the only simple geometrical form of 
finite array with the simple AFM structure which can 
be bordered by dense lines for particles with alternating 
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magnetic moments. Here the three possibilities of 
orientation of rhombus with respect to the lattice 
corresponds with three alternating AFM structures 
present at Fig. [TJ For this reason, at zero magnetic 
field only the rhombus could provides ideal AFM 
ordering. Numerical analysis confirm this suggestion 
demonstrating simple two-sublattice structure, with lines 
from one sublattice inclined by the angle 60° to the 
border lines, for rhombuses with different sizes and 
shapes, see an example on the Fig. [5] 



«2-4» rule is broken for the central site of hexagonal 
sample. In the nearest-neighbor approximation, the 
effective field on the central site equals to zero, and 
the central magnetic moment is very sensitive to the 
application of the magnetic field. For large hexagonal 
sample with real magnetic dipole interaction, the field 
of the turn of this moment is also small, 

H hex = -( 6( r 6 V m +3a+3e-i)^ = 0.428266^ , (21) 
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Phc. 9: Minimal configuration for a rhombus-shaped finite 
system at zero magnetic field. 



and the magnetization reversal is starting by the turn 
of the central magnetic moment. But detailed analysis 
shows that the field for the turn of the moments 
neighboring the center is growing fast with the distance 
to the center, and the magnetization reversal is following 
the same scenario as for infinite or semi-infinite systems, 
see the next subsection. 



2. triangle and hexagon 



3. rectangle 



The samples shaped as a right triangle or a hexagon 
with ideal AFM structure cannot satisfy the above 
condition, an absence of the border line from one 
sublattice. This leads to appearance of some kind of 
domain structures in the ground state, see Fig. [TU1 

For triangular sample, the domain structure with two 
AFM domains turned at the angle of 60° , is present in the 
ground state. For the hexagon, the minimal configuration 
consists on three domains, where the sublattice lines are 
inclined on 120°. The three domain walls in this structure 
are going out from the corners and coming together at 
the center of the sample in such way that the structure 
has the three-fold symmetry axis, see right part of the 

Fig.nni 

It is worth noting here that the corresponding 
domain walls have zero energy in the nearest-neighbor 
approximation, but the energy losses caused by the walls 
are small for long-ranged dipole interaction as well. For 
these samples most of the «bulk» sites have the optimal 
«2-4» neighborhood, with the only one exception: The 
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Phc. 10: minimal configurations for the triangle with the side 
21 (left) and for the hexagon with the side 8 (right). 



For trivial geometrical conditions, the rectangular 
sample can have only one pair of sides (two parallel sides) 
parallel to elementary translation vector; these sides have 
standard AFM structure with interlacing of up and down 
magnetic moments, see Fig. 1111 Two other pair of borders 
of rectangular are parallel lines of sites with relatively 
low density within the line (the distance between the 
sites equals to a-\/3 ~ 1.73a), but with a small enough 
distance (a/2)^/3 < a with the next equivalent lines. The 
preferable type of ordering for two such lines corresponds 
to the parallel alignment of the magnetic moments within 
the line with the antiparallel orientation for neighboring 
lines («2-0» type surface state) whereas the ideal AFM 
structure leads to much less favorable «1-1» type surface 
state. For this reason, one can expect the deformation 
of the ideal AFM structure in a rectangular sample. The 
numerical analysis confirms this simple speculation based 
on nearest-neighbor approximation, see Fig. 111! 

The ground state of rectangular sample contains 
zigzag deformation of the characteristic lines of the 
ideal AFM structure such that the ideal AFM sign- 
interlacing distribution is kept at the dense border, 
with the specific structure of finite half-hexagons on 
the less dense borders. This structure provides the 
aforementioned picture of saturated lines parallel to less 
dense borders with «2-0» type surface state. The zigzag 
lines are always parallel to the less dense borders, whereas 
the hexagon surface structure depends on how the long 
side of rectangular sample is oriented respectively to 
the elementary translation vector, see Fig 



square sample. 



11 a 



Fig. 11(c) for elongated rectangles and Fig. 11(b) 



and 
for a 
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Phc. 11: Minimal configurations for rectangular samples 
of different shapes; (a) — rectangle, elongated along the 
elementary translation vector; (b) — square sample, (c) 
- rectangle, elongated perpendicularly to the elementary 
translation vector. 




Phc. 13: (color online) Magnetization function for the finite 
systems; hexagon with the side length 8 (215 sites) and square 
12 x 12 (144 sites), in comparison with that for infinite system. 
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Phc. 12: Ground state magnetic configuration for the circular 
sample of the radius 9 found by Monte-Carlo simulations. 



4- circle 

The circular samples are considered here as an 
example of systems with non-small surface roughnesses. 
The numerical analysis demonstrates the essential 
deformations of the ideal AFM structure for such 
samples, see Fig. [12] Quite non-regular border structure 
having sites with 3 or 4 nearest neighbors, leads to big 
variety of the local surface configurations with different 
numbers of parallel and antiparallel nearest magnetic 
moments, like «2-l», «3-0» and «3-l» local states. But 
note that the local AFM structure, in particular, «4-2» 
condition for the local states, is kept for most bulk sites. 



B. Magnetization processes for finite samples. 

The magnetization functions for finite systems: 
hexagon with the side length 8 (215 sites) and square 
12 x 12 (144 sites) are presented on the Fig. [T31 The 
result for infinite system is cited for comparison. 

The magnetization function for the finite arrays in 
many details repeats such a function for infinite system. 
All typical fields of the structure change are reduced due 
to long-range nature of the dipole-dipole interaction and 



its cropping on the length of the order of the sample 
size R. This feature is more strikingly expressed for 
the saturation field, where the respective sum converges 
with the distance as R~ l , rather than for the field of 
antiferromagnetic destruction, which evaluates by the 
rapid convergent alternating sums. At small fields the 
same the topological domain wall scenario is present, 
with the orientation of the domain walls consistent with 



the shape of the sample, see Fig. 14(a) Step-like behavior 



(plateau behavior) of the function of the infinite system 
corresponding to the becomes less pronounced in the 
case of the finite array due to the incommensurateness 
between the superlattice structure and the shape of the 
sample. However the period is small for the triangular 
superlattices of the small parameters (V3 h 2) and 
incommensurateness does not have impact yet. In this 
case typical plateaus are visible on the magnetization 
curve and additional fine structure of the form of the 
low steps appears on them. This fine structure is relating 
to the change of superstructure contours but not to its 
reorganization. 

As has been shown in Sec. [Vj the edges of finite array 
are "weak points" for a single-dot instability, and the 
values of the instability field for a dot at the vertex 
of edges are lower than in the bulk, see Eq. (JTSJ) and 
Eq. (|19p . Thus for a finite sample with the increasing 
of the external field the saturation starts at the border 
at lower fields and the area of the superlattic e will b e 
reduced. Such behavior is shown on the Fig. |14(b)| h 
Fig. 14(c) As one can see from the Fig. [HI configurations 
for the hexagon are in the good correspondence with the 
configuration for the infinite system. The only difference 
is the presence of saturated dense border lines at high 
enough fields. 
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(b) h = 2.0 



(c) h = 6.0 
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(d) h = 6.1 



Phc. 14: Magnetic structure for hexagonal sample with the 
size 7 at some characteristic values of the magnetic field. 



VII. DISCUSSION AND CONCLUSIONS 

Let us now discuss the general regularities revealed in 
the behavior two two-dimensional systems with dipolar 
coupling of Ising magnetic moments, the triangular 
lattice and the square lattice of magnetic particles. First 
note that for nearest-neighbors interaction the features of 
ordering for these two lattices are of principal difference; 
the square lattice Ising moments at finite temperature 
shows typical phase transition to the phase with the long 
range chessboard AFM order, whereas the AFM Ising 
triangular lattice system still unordered until T — > O i 33 i 34 
The reason, well discussed in the literature, is based on 
the frustrated character of triangular magnets with AFM 
interaction.— The nearest-neighbors Ising AFM system 
is an extreme example of such feature, because it has 
topological defects with zero energy, which can lead to 
destruction of long-range order i 35 ! 36 

For our case of long-ranged dipolar interaction the role 
of frustration is not so crucial; it is enough to mention 
that for a square lattice of Ising moments with dipolar 
coupling the next-nearest-neighbors in AFM structure 
are frustrated as well. On the other hand, the topological 
defect lines having zero energy for triangular lattice 
nearest-neighbor system gain finite (but small) energy 
for real dipolar interaction. But the effects of frustration 
for nearest-neighbors interaction produce an essential 
difference in the behavior of the systems of interest, 
magnetic dot arrays with square and triangular lattices. 

First let us mention the properties common for both 
square lattice and triangular lattice. Their ground state 
at zero field is antiferromagnetic, and the saturated 
(ferromagnetic) state is present for high enough magnetic 



fields. For high magnetic fields, nearly the saturation, 
the behavior of both lattices is practically the same. 
The destruction of the saturated state for both lattices 
happens through creation of a superlattice of reversed 
magnetic moments, and at h — > h sa t both superlattices 
are triangular. We can only mention a quantitative 
difference of the saturation fields, the value for triangular 
lattice, 11.034too/o 3 , is a bit higher than for a square 
lattice, -ffsat, square = 9.033622mo/a 3 . This difference just 
reflects the fact that the square lattice is less dense than 
triangular. But for the region of low fields, corresponding 
to the destruction of AFM state, and for intermediate 
fields, where non-saturated states with (to) = (0.2- 
0.8)mo occurs, the behavior of these two lattices is 
completely different. 

For the square lattice the states with small, but non- 
zero (to) can be obtained from the chessboard AFM 
structure by reversing the magnetic moments of a small 
fraction of magnetic dots upward, leaving the remainder 
undisturbed*^ For a triangular lattice, the destruction 
of the ideal AFM structure is determined by creation 
of the system of topological lines (domain walls). The 
creation of such a single line (or small density of such 
lines) effects approximately one half of the particles in 
the system. Thus, this phenomena appears at much lower 
field (compare the values 2.646TOo/a 3 and 0.7086too/o 3 
for these two lattices) but needs to overcome much higher 
potential barrier. 

For intermediate region of fields, the common property 
for both lattices is that the ground states are mostly 
characterized by complicated superlattices with different 
densities of "up" and "down" magnetic moments. But 
the process of magnetization for an array with square 
dot lattice is influenced mainly by an interplay between 
square and (distorted) triangular superlattices of "down" 
moments. Such behavior is dictated by a compromise 
between the optimum of interaction energy (evident for 
triangular superlattice) and the absence of the distortion 
energy for square superlattice commensurate with the 
underlying lattice. In contrast, for a triangular dot 
array triangular superlattices are optimal from any point 
of view. Such structures occur at the main regions 
of field, producing well-defined "shelves". However, the 
frustration effects are present for any triangular lattices 
and superlattices, leading to complicated character of 
the transition between the states with superlattices of 
different period, especially small periods like ay/3 and 
2a. 

It is also reasonable to discuss briefly problems that 
remain beyond the scope of this paper. Of course, the 
problem of the transition from the out-of-plane (Ising) 
states to planar states when the magnetic field varies 
for finite anisotropy for a single dot is of interest. 
Anisotropy energy for n-th dot can be written as 
(1 /2)moH lin (a n ■ el) 2 , where the anisotropy field, H an 
characterizing the strength of the anisotropy energy, is 
introduced. In fact, this means the construction of a 
phase diagram on the H, H an plane. One can expect that 
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Ising states are stable for sufficiently large anisotropy, 
when the anisotropy energy is comparable with the 
energy of interaction between neighboring dots, H an > 
mo /a 3 . The analysis of the square lattice corroborates 
this rule^ and demonstrates that the Ising states 
are stable at H an > 5mo/a 3 . It is rather difficult 
to analyze planar states analytically. Our preliminary 
numerical data indicate that the non-saturated states 
at small anisotropy correspond to complex noncollinear 
structures, which are characterized both by significant 
two-dimensional inhomogeneity with a scale of about 
the sample size, and by the presence of regions 
where neighboring magnetic moments are substantially 
noncollinear, and full description of such states can be 



done only numerically. The complete analysis of non- 
collinear states is evidently far beyond the scope of our 
article. 
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Appendix. Numerical method description 

To find the global minimum of the energy of the system 
the Monte-Carlo approach with a simulated annea ling 
(MC-SA) method was used, see original articles Refs. [40, 
SU and the textbook Ref. |H. For MC-SA realization, 
first, the random initial configuration was selected. Every 
iteration of the MC-SA method consists of N moment 
reversal attempts on the random site, where N is the 
site number in the sample. The main idea of simulated 
annealing is that the probability of the reversal is non- 
zero even if the energy is growing after this reversal; 
otherwise, the system with a high probability will be 
"frozen" in some local minimum. The probability depends 
not only on energy gain for reversal, which equals to 
Hm (H is a value of the field created by other dots 
on a given ) but also on a global time-varying parameter 
T called the temperature. If the reversal is favorable in 
energy, the moment is always reversed, irrespective of the 
temperature. But even if the reversal is unfavorable, the 
non-zero probability of reversal is chosen as follows: flip- 
over takes place if Htuq < T\ logp|, where T is the current 
value of temperature, p is a random value generated in 
the range < p < 1. Here the parameter temperature 
determines the strategy of the minimization: for large 
T, the evolution is sensitive to coarser energy variations, 
while it is sensitive to finer energy variations when T is 
small. Thus the meaning of the temperature is the same 
as for annealing in metallurgy involving initial heating 
and controlled cooling of a material thereby avoiding 
defect formation. 

The temperature is changing according to the quantity 
of full steps of MC-SA of the sample n as follows, 
T = Tq min[«;, (no/n) a ]. Here the parameter Tq was 



chosen as (0.2-0.4)toq/ci 3 , and the cutoff parameter k was 
equal to k — 3.0 such that the initial temperature was 
high enough compared with the interaction energy. The 
optimal values of other parameter are defined by the trial 
runs; no was equal to 10 4 - 5 • 10 4 and the value of the 
a-index was taken 1/4 or 1/5. 

The temperature decreases with the process evolution 
and magnetization reversals take place more rarely 
with decreasing energy. The process was stopped if 
the energy did not become less than the previous 
minimum during the previous lOrio iterations. Then 
the next random configuration was generated and the 
process was repeated. Looking at the energies of minimal 
configurations for every run one can estimate the 
probability that such a configuration corresponds to 
the true minimum. For example, almost all processes 
of the cooling go through the same minimum far from 
the regions of the parameters which correspond to the 
transition between different states. 

For the configurations on the infinite lattice we took 
the rhombuses with the consecutive increasing periods up 
to 16 x 16 and used periodic boundary conditions. The 
energy of the configuration was recalculated to the energy 
per one dot thereby making possible the comparison 
of the results for the different periods. A rhombus was 
selected to admit the maximum possible configuration set 
to consideration, but minimal configurations often have 
higher symmetry. 

A similar approach was used to build a magnetization 
function. Firstly, all ranges of the fields were passed with 
large intervals and small periods of the rhombuses. After 
that the regions of the structure rebuilding were detected, 
where the field intervals were taken smaller and periods 
were taken larger. 



